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Abstract 

A linearization procedure is proposed for Ermakov systems with 
frequency depending on dynamic variables. The procedure applies to 
a wide class of generalized Ermakov systems which are linearizable in 
a manner similar to that applicable to usual Ermakov systems. The 
Kepler-Ermakov systems belong into this category but others, more 
generic, systems are also included. 


1 Introduction 

Ermakov systems have merited special attention in recent years. A 

trend in the latest developments on the subject, is to focus attention in 
some special features of subclasses of Ermakov systems. These subclasses 
are frequently more flexible and may be tailored to suite some particular 
application or special purpose. Among others, we quote applications such as 
the identiheation of the Hamiltonian character in special circumstances iJ; 
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the determination of a second constant of motion in certain other particular 
cases PI; the identihcation of the structure of associated Lie group of point 
symmetries 1,0 and the extension of the Ermakov systems concept itself 
to higher dimensions ||T0|-||T3 


A central feature of Ermakov systems is their property of always possess¬ 
ing a hrst integral, commonly known as the Lewis-Ray-Reid invariant. The 
Lewis-Ray-Reid invariant can be used to construct nonlinear superposition 
laws, another important feature of Ermakov systems concerning their general 
solutions |]T4 . 


Recently, Athorne shown that Ermakov systems in their usual form, where 
the frequency function depends only on the time variable, are linearizable 
r5|. In the linearization process, the Lewis-Ray-Reid invariant plays a cen¬ 


tral role, analogous to that played by the angular momentum in the lineariza¬ 
tion of the two-body problem of classical mechanics |]^. The linearization 
process, besides being an interesting mathematical attribute in itself, has 
found some important applications. The linearized form of the Ermakov sys¬ 
tems plays a central role, for example, in the resolution of a problem arising in 
the theory of two-layer, shallow water waves [0, in the singularity analysis 
of Ermakov systems [0, 0 and in the study of higher dimensional Ermakov 
In such applications, however, a restricting condition on the 
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systems [| 

class of Ermakov system is that its frequency function depends only on time, 
a fact that imposes limitation on the scope of the linearization process. In 
a recent study, Athorne has found a class of dynamical systems, the so 
called Kepler-Ermakov systems, that generalizes the usual Ermakov systems 
while preserving the property of being amenable to linearization. The natu¬ 
ral question to ask at this point is, therefore, whether the Kepler-Ermakov 
systems are the only linearizable generalization of the Ermakov systems or 
if there exist other perturbations of usual Ermakov systems that are also 
amenable to linearization. 

The frequency function entering the Ermakov system may depend on 
the dynamic variables and their derivatives, without any restriction on the 
existence of the Lewis-Ray-Reid invariant [0|. The resulting generalized 
Ermakov systems, with frequency function depending on the dynamic vari¬ 
ables and their derivatives, besides time, have appeared in the study of a 
generalized time-dependent sine-Gordon equation 
Hamiltonian character |^, the existence of a second exact invariant ^ 
the Lie point symmetry group |p of Ermakov systems have been analyzed 
for such generalized frequency functions. In the present work we address 


More recently, the 
and 
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the question raised at the end of the previous paragraph and study the lin- 
earizability of generalized Ermakov systems where the frequency function 
may depend on the dynamic variables and their Erst derivatives, besides 
time. Our motivation is to hnd a more general class of dynamical systems 
still amenable to linearization, on the same line as the usual Ermakov or 
Kepler-Ermakov systems. In fact, as shown in detail in section III, the 
Kepler-Ermakov dynamical systems are nothing but a special class of the 
linearizable generalized Ermakov system. This result puts the linearizabil- 
ity property of Kepler-Ermakov systems in a novel, more generic and sound 
basis. It must be stressed, however, that the class of generalized linearizable 
Ermakov systems is much ampler then the class of Kepler-Ermakov systems, 
as will be seen in section II. 

The paper is organized as follows. In section II, a linearization proce¬ 
dure is proposed for generalized Ermakov systems. The class of frequency 
functions for which the linearization procedure applies is obtained in terms 
of a wide category of generalized Ermakov systems that can be linearized. 
In the continuation, it is shown that the general solution for the nonlinear 
equations can be recovered from the solution for the corresponding linearized 
equation. The linearization of usual Ermakov systems is recovered as a spe¬ 
cial case of the more general theory. In section III, Kepler-Ermakov systems 
are shown to belong to the class of generalized linearizable Ermakov systems. 
The linearization procedure for the Kepler-Ermakov system is illustrating in 
the case of a particular non-central force problem. In section IV an example 
of a generalized linearizable Ermakov system is presented which is not of 
the Kepler-Ermakov type. The linearization of this system is also exhibited. 
Section V is dedicated to the conclusions. 


2 Linearization of generalized Ermakov sys¬ 
tems 


in 2-D, is a pair of coupled, nonlinear second 

order differential equations. 


An Ermakov system |M 14 


x + u'^x = -^f{y/x), 
yx^ 

( 1 ) 

y + uj‘^y = -^g{x/y), 
xy^ 

( 2 ) 
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where / and g are arbitrary functions of their indicated arguments and a;, the 
so-called frequency function, may depend on time, the dynamic variables and 
their derivatives [Q. In other words, oj = a;(t, x, y, i, y, x,y,...). Either for 
physical or else for simplifying reasons, we restrict considerations to the cases 
where the frequency function depends at most on the velocity components. In 
this work, the nomenclature “usual Ermakov system” is used when a; = a;(t), 
and “generalized Ermakov system”, or Ermakov system for short, is used 
when the frequency has a more general dependence. 

Independently of the way on which a; depends on the dynamic variables, 
the system (|I|-^) always possess the Lewis-Ray-Reid invariant 

^ = ^ixy- yxf + U{y/x ), (3) 


where U is dehned by 

U(vlx) = f(X)dX + f'’ g(X)dX. (4) 

The invariance of I can be directly verihed by checking that dl/dt = 0 along 
any trajectory of the Ermakov system. 

For the specihc purpose of this work, polar coordinates x = r cos 9 and 
y = r sin 6 are more appropriate. In this system of coordinates the Ermakov 
system becomes 


f — + = F{6)/r^ , (5) 

re + 2x9 = -{dV/de)/r^, (6) 

where V{9) = U{ta.n9), F being dehned in terms of / and g by 

f{tan9) + g{cos9) 

In polar coordinates, the Lewis-Ray-Reid invariant becomes 

/ = + V(e). (8) 

As already pointed out before 0a. the concept of a generalized Ermakov 
system has originated from the observation that uj may depend arbitrarily 
on the dynamic variables and that, as a consequence, only two and not three 
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arbitrary functions are necessary to specify the Ermakov systems. Indeed, 
redefining 

o;^ I—o;^ — F(6')/r^ (9) 


is equivalent to absorbing F into cn, that is, to making F = 0 in the pair of 
equations (||-||). Thus, strictly speaking generalized Ermakov system con¬ 
tains only two and not three arbitrary functions. This redefinition, however, 
is not essential in the present context and the usual notation is adopted in 
order to easy the comparison between our results and the results found in 
the literature. 

The linearization procedure for usual Ermakov systems is described in 

TB. 


ref. 


Our strategy for the linearization of generalized Ermakov systems 


follows the same spirit of the procedure proposed for usual Ermakov sys¬ 
tems. However, the assumption of a more generic dependence in uj suggests 
a straightforward generalization of the techniqne which leads to interesting 
new results. In this way, we are able to hnd classes of generalized Ermakov 
systems, specihed by appropriate dependence of uj in the dynamic variables, 
that are linearizable. For this purpose we follow Athorne [|T^ and introdnce 
the new dependent variable V’ dehned by 


V’ = p(t)/r, (10) 

where p{t) is an arbitrary, unspecihed fnnction of time. The new independent 
variable will be chosen as the angle 9. In order to obtain the equation of 
motion in the new variables, it is necessary to express the time variation of 
9 in terms of the invariant. This is achieved by use of the Lewis-Ray-Reid 
invariant (§), from which we construct 

9 = h{9-I)/r\ (11) 


where a function h, parametrically dependent on the nnmerical valne of /, 
was dehned by the relation 

h{9;I) = V2iI -V{9)Y/^ . (12) 

At this point it is interesting to compare eq. ([Tl|) with the corresponding 
eqnation for 9 in the two body problem. Here, the Lewis-Ray-Reid invariant 
plays the role of angular momentum whereas V{9) is a specihc feature of the 
non-central natnre of the motion. 
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Now using the Ermakov system (§-D and equations we easily 

arrive at the transformed equation 


h^ej)^+h(ej) 


dh{9] I) dip 

86 le 


f {hHe; n + F(B)) i, 


pHp + 

Ip^ 


(13) 


In order to obtain a linear equation in transformed variables, the dependence 
of u in the dynamic variables must be chosen such that the right hand side 
of eq. (^) becomes a linear function of ip and dip/dO. This condition yields 
the compatibility condition for linearization, 




(14) 


where a, b and c are arbitrary functions of the indicated arguments. A 
dependence on the numerical value of the Lewis-Ray-Reid invariant I was 
included for maximal generality. Notice how the extra dependence in uj allows 
for a more general solution. In the case of a; = uj{t) treated by Athorne, only 
p + = 0 with a = b = c = 0 was consistent with linearity. In the 

present case, p remains arbitrary and three new functions a, b and c are 
introduced, leading to more general expressions for uj which are compatible 
with linearization. The resulting frequencies are given by 


UJ 


= -'^ + 

p 


dip 


ip^ 


a{e-,i)-^ + b{e-,i)ip + c{e-,i)]^. 


dd 




(15) 


In principle, we could include higher order derivatives in the right hand side 
of equation (|l^ , without restrictions to linearity. This possibility, however, 
was neglected, mainly for simplifying reasons. 

Now restricting the frequency function of the form (|I^), we hnd the lin¬ 
earized Ermakov equation 


h\e-j) 


d^ip 


dh 


dip 




dO 


{h\e-I) + F{9)-b{9-,I))iP = c{9-I). 


(16) 


Frequencies not of the form (P!5|), necessarily correspond to Ermakov equa¬ 
tions that are not linear when the transformation (r, t) —>■ {ip, 9) is performed. 

In order to express the frequencies of the linearizable systems in terms of 
the original polar coordinates, we can use 


dilj/d9 = ip/9 = —{pr — pr)/r‘^9 , 


(17) 
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so that 


w = -- + 
P 


p {pr — pr) 


pr^ 


A{9,r^9) + 


B{9,r^9) C(9,r‘^9) 


+ 


pr'- 


where 


A(9,r^9) = -a{9]I)/{r^9), 

B(e,rH) = *(»;/). 

C(e,r‘e) = c(9;/) 


(18) 


(19) 

( 20 ) 
( 21 ) 


are new arbitrary functions. Notice that the dependence of A, B and C in 
equations is of the correct type, since the invariant (|]) depends only 

on 9 and r‘^9. 

The resulting linearizable generalized Ermakov system in polar coordi¬ 
nates has the form 


r 


r9 -|- 2r9 


r9^ -^r+ f^^ A{9, r^9) + 

p pr^ 

B{9,r^9) C(9,r^9) _ F(9) 

r3 pr'^ ’ 

dVld9 


( 22 ) 

(23) 


and involves six arbitrary functions, namely p, A, B, C, F and V. Other 
classes of generalized Ermakov systems amenable to linearization can exist, 
besides those generated by a dependence of uj on higher order derivatives. We 
may convince ourselves of this by noticing that ("0, 9) is only one particular 
choice of transformed variables. Other different choices can be used, possibly 
leading to other classes of linearizable generalized Ermakov systems. The 
great generality of the system (^^23), however, is already sufficient for the 
purposes of this paper. 

Some additional remarks are in order here. First, as already pointed out, 
the function F can be dropped from the system by an adequate redehnition 
of oj. This is manifest in equation (|22|) where we see, by inspection, that 
B can account for F. The second and perhaps more interesting remark is 
that the linearizable generalized Ermakov system ([2^-^3D can be expressed 
in autonomous form, by means of a quasi-invariance |^, ^ transformation 


r = r / p , 9 = 9 



(24) 
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( 25 ) 


In these variables, we have the autonomous representation 


r9'^ + — {Ar'+ B/r + C) = 


F{9) 


f9" + 2f'9' = 


dV{9)/d9 

^3 


(26) 


where the prime stands for derivative with respect to i and A, B and C 
are functions of f‘^9' and 9. The time dependence in (p2|-p3|), is thus, in 
a sense, spurious. Consequently, there remains, in fact, only four funda¬ 
mental arbitrary functions in the linearizable generalized Ermakov system 
([2^-p^, that is, F may be incorporated in C and p may be eliminated by 
the quasi-invariance transformation. It must be stressed, however, that the 
quasi-invariance transformation is not an essential step in the linearization 
procedure. 

We may now recover the solution for the nonlinear dynamics from the 
solution for the linearized Ermakov system. Let 


-0 = ip{9] /, Cl, ca) = ciil)i{9) + C2ip2{9) -f il)p{9) 


(27) 


be the general solution for the linearized equation (|^), where Ci and Ca are 
arbitrary constants, 'ipi and '0a are two linearly independent solutions for 
the homogeneous part of ([T^), and 0p is any particular solution for (pT]). 
The solution for the nonlinear system follows from the dehnition (p!0|) of the 
linearizing variable 0 and relation (^), from which we hnd 

9 = h{9■I)^p\9■I,c^,C2)/p\t). (28) 


This is a separable hrst order, ordinary differential equation, equivalent to 
the quadrature 

dX A dX 

J h(A;/)02(A;/,ci,ca) J p‘^{X) 

where J is the fourth integration constant of the equations of motion. The 
quadrature (|^) locally yields the equivalent relation, 

9 = 9{t;I,J,Ci,C2), (30) 

which is the general solution for the angular variable, involving four arbitrary 
integration constants, or implicitly 



t = t{9-,I, J, Ci,C2) . 
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( 31 ) 








To obtain the orbits, we use the dehnition of ip, the solution m for the 
linearized system and the relation (^T]), 


r = r(6';/, J,Ci,C 2 ) = 


p(t(0; J, J,Ci,C2)) 


(32) 


i/>(6';/,ci,C2) 

The time evolution of the radial variable follows from the dehnition of ip, 
the solution for the linearized Ermakov system and the relation (|30|), from 
which we obtain 


r = r(t; J, J, Ci,C 2 ) 


pit) 

tpit;I, J, Ci,C2) ’ 


(33) 


where 'ip(t]I,Ci,C 2 ) = ip {6[t] 1, 01 , 02 )] I 02 )- While equation (|32|) repre¬ 
sents the orbits for the nonlinear generalized Ermakov system, equations 


and (^) give the time evolution of the dynamic variables, in terms of 
the general solution for the linearized system. 

A particular case in the class of generalized linearizable Ermakov systems 
(|2^-^^, are the usual Ermakov systems, with frequency depending only on 
time, uj = ui{t). Let A = B = C = Q in the dehnition ( p!8D of generalized 
frequencies, which implies 


p -f a;^(t)p = 0 . 


(34) 


Equation (^) is the equation for a harmonic oscillator p, P| with time- 
dependent frequency. In this particular case of usual Ermakov system, the 
linearization ([I^) reduces to the form 


I'pSi I)^ + h{0; I) 


dh{9; I) dip 

de He 


+ (h^(D;!) + F(D)) ^. = 0, 


(36) 


which is a homogeneous linear equation. This result is in full agreement with 
that of At home II , although in a slightly diherent form, due to our par¬ 
ticular choice of linearizing variables. Illustrative examples of linearization 
of usual Ermakov systems can be found in |T^, In the following 

section, a non-trivial generalized Ermakov system that was not, until now, 
treated as a member of the class of Ermakov systems is examined and its 
linearization explicitly calculated. 
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3 Kepler—Ermakov systems 


Kepler-Ermakov systems [20| are given by 


= F(9)/r=‘-G(9)/^^ 
rB + 2re = -(dV(e)/d0)/r^, 


(36) 

(37) 


with F, G and V are arbitrary functions of the indicated arguments. Kepler- 
Ermakov systems were introduced as a linearizable perturbation of usual Er¬ 
makov systems. The term G in equation (|36D destroys the usual Ermakov 
character of the system. However, the system (^6|-^) belongs to the class 
of linearizable generalized Ermakov systems derived in section II and con¬ 
sequently, qualifies as a generalized Ermakov system that is linearized by 
the procedure outlined above. This result sheds new light on the cause why 
Kepler-Ermakov systems are linearizable. In fact, choose 


A = B = 0 , G = G{9) , p=l (38) 


in the definition of linearizable generalized Ermakov system (p^p3|). This 
choice puts the system in the form of a Kepler-Ermakov system and the 
corresponding linearization is given by 


which, unlike the linearization of usual Ermakov systems, is a non-homogeneous 
linear equation. The non-usual character of the frequency function associated 
to Kepler-Ermakov systems is manifest in the equation 


UJ 


2 


G{9)/r^ 


(40) 


which, evidently, shows a dependence on the dynamic variables. 

A typical example of linearization of a Kepler-Ermakov system can be 
found in reference . Other relevant systems that can be cast in the Kepler- 
Ermakov form are some of the superintegrable systems treated by Winternitz 
et al. 1^ and by Rahada |^. To inspect one of these systems as a Kepler- 


Ermakov system and to show how the linearization process works in a specihc 
case, let us consider the Hamiltonian 


H 



/fo 1 

rp ' p2 


/ Qi + g 2 cos 9 
\ sin^ 9 



(41) 
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where /xq) 9i) 92 and are positive constants. When = 0, the correspond¬ 
ing Hamilton-Jacobi eqnation is separable in parabolic and polar coordinates, 
a fact that results from the underlying dynamic symmetry algebra |^. For 


g 3 7 ^ 0, however, the system is still completely integrable, as shown in 
Straightforward computation of the canonical equations 


r = dH/dpr , 0 = dH/dpe , pr = —dH/dr , pe = —dH/dpe 

shows that the Hamiltonian (^T]) yield a Kepler-Ermakov system, with 

V{6) = {gi + g2Cos6)/siii^ 6 , 

F{e) = 2{V{e)+g^), 

G{e) = /io 

in the equations of motion (|3^37). The corresponding linear equation 
is easier to treat in terms of a new time parameter T{6] /) dehned by 


T{9;I,J)= I 


0 dX 


h(A;/) 


+ J • 


(42) 


(43) 

(44) 

(45) 


(46) 


where J is an arbitrary constant. Interestingly, this new time parameter de¬ 
pends parametrically on the numerical value of the Lewis-Ray-Reid invariant 
J. In terms of the new independent variable, the linearization becomes 

+ {I + 93)i’ = ho , (47) 

the equation for a harmonic oscillator with a time-independent driving, 
which is always exactly solvable. The explicit form of the solution will de¬ 
pend on the parameters I and g^. To construct a particular case, take / > 0 
and 5^3 > 0. This choice together with (^71) yields, 

= ^jJ{T; Cl, C 2 ) = Cl cos(/ gsY^'^T -f C 2 sin(/ -h gsY^'^T , (48) 


where ci and C 2 are arbitrary constants. To express V’ in terms of the angle 
6, we use formula (^) which, for / > 0, furnishes 


j) = 



2 / cos 9 + g 2 

(j| + 4/(/-<,i))i/2_ 


+ J . 


(49) 
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After substituting this result in (^), we find ip = ■^(6*; I, J, Ci, C2). The orbits 
are then given by 

r = r(6'; /, J, Ci, Ca) = I, Ci, C 2 ) , (50) 

which is obtained from the dehnition of the linearizing variable pj and the 
fact that p = 1 in the case of Kepler-Ermakov systems. The actual form 
of the orbit is intricate and will be omitted here. In spite of the orbits 
being analytically given in the present example, the time evolution of the 
dynamic variables cannot be given analytically because the quadrature (pOD 
is not expressed in terms of elementary functions. Similar reasoning apply 
for other values of the parameters I and g^. 


4 Linearization and reduction to free motion 


In this section, we reverse the arguments and search for classes of Ermakov 
systems whose linearized form is simply a free motion. Among those, we 
hnd examples of linearizable Ermakov systems that are not of the Kepler- 
Ermakov type. In addition we hnd that a whole class of Ermakov systems, 
still depending on two arbitrary functions, are reducible, by the linearization 
process, to simple free motion. Let 

a = hdh/d0 , b = + F , c = 0 (51) 

in the linear equation ( 0 ). yielding free particle motion, 

= 0 . (52) 

In this case the general solution is 

p) = Cl+ C 26 , (53) 


where Ci and C 2 are arbitrary constants. In this case the quadrature (^) has 
the form 

J h(A;/)(ci + C2A)2 J p2(A) - • ^ ^ 

For any given function h{ 6 -, I), we locally hnd, from equation (p^ , that either 
t = t{6] /, J, Cl, C 2 ) OT 6 = 6{t] J, J, Cl, C 2 ). The corresponding orbits are given 
by 

p(t(d;/, J,Ci,C2)) 


r = r(6';/, J,Ci,C 2 ) = 


Cl + C29 


(55) 
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whereas the time evolution of the radial variable follows from 


r = r(f; J, J, Ci,C 2 ) = 


p{t) 


Cl + C 26 '(f;/, J, Ci,C2) 


( 56 ) 


This completes, in the quadrature sense, the integration of the equations of 
motion. 

Performing the steps in reversed order, we hnd, from (pT H^) and 
the dynamical system in the original variables 


A = {dv/de)/{r% , B = {r^ey + F{e) 


C = Q 


(57) 


This yields the frequency function 


2 p ^ {pr — pr)dV/d9 ^ + F 

p pr^e 

and the corresponding linearizable generalized Ermakov system 


(58) 


pr — 'pr 
rO + 2r9 


{pr — pr) dV{9) 

r‘^9 d9 
dV{9)/d9 


(59) 

(60) 


As a hnal step, we can express the dynamical equations in cartesian form 
dm), with 


UJ 


2 


P ( xy-yx \ 
p \x^ + 1)“^ j 

9{x/y) = 


/ {px - px)x + {py - py)y 
y px‘^y‘^{xy — yx) 
-f{y/x). 


f{y/x), (61) 
(62) 


These systems which are of the generalized Ermakov type certainly have 
remarkable properties: despite having two arbitrary functions, p and /, their 
solution can be reduced to the equation of a free particle plus the quadrature 

0 )- 


5 Conclusion 

In this paper, an extensive class of generalized Ermakov systems, character¬ 
ized by several arbitrary functions was shown to possess equivalent linear 
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form. The linearization process depends basically on the existence of the 
Lewis-Ray-Reid invariant, a fact that is always assured for Ermakov sys¬ 
tems, both in its usual or generalized form. As a particular result of the 
theory, it was shown that the linearization of Kepler-Ermakov systems is 
nothing but a special case of the more general linearization of some classes 
of generalized Ermakov systems. While usual Ermakov systems are always 
linearizable, the same does not apply to generalized Ermakov systems. In the 
present work, use was made of the linearization variables {'4^,0), the choice 
that originates the system (^^- |^) - Different choices of linearizing variables 
could eventually lead to other classes of linearizable generalized Ermakov 
systems. As a hnal comments we remind that linearization could be a use¬ 
ful criteria to qualify or classify Ermakov systems. For instance, using the 
linearization, Athorne has classihed usual Ermakov systems with respect to 


the rational character of their solutions 


The same strategy can be 


applied to linearizable generalized Ermakov systems. This could be adopted 
as an alternative to the symmetry criteria proposed by Leach [10| . At least 


the class of linearizable Ermakov systems seems more extensive - four arbi¬ 
trary functions - than the class of generalized Ermakov system with Lie point 
symmetry which is specihed in terms of only two arbitrary functions. 
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